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Abstract 

The green microscopic alga Ostreococcus tauri has recently emerged as 
a promising model for understanding how circadian clocks, which drive the 
daily biological rythms of many organisms, synchronize to the day/night cy- 
cle in changing weather and seasons. Here, we analyze translational reporter 
time series data of its central clock genes CCA1 and TOC1 for a wide range of 
daylight durations (photoperiods). The variation of temporal profiles with day 
duration is complex, with the two actors tracking different moments of the day. 
Nevertheless, all profiles are accurately reproduced by a simple two-gene tran- 
scriptional loop model whose parameters are affected by light only through the 
photoperiod value. We show that this non-intuitive behavior allows the circa- 
dian clock to combine flexibility and robustness to daylight fluctuations. 

Keywords: circadian clocks; robustness; flexibility; Ostreococcus tauri; 
entrainment 
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1 Introduction 



Most biological functions are controlled by complex networks of molecular inter- 
actions which routinely achieve sophisticated information processing and decision 
making. These networks, featuring manifold feedback and feedforward regulations, 
display a highly nonlinear collective dynamics Q3-[3]|. A natural idea is that besides 
reproducing observations and checking the consistency of a biological hypothesis, 
mathematical modeling can help us to identify design principles at work in cellular 
machinery flU. For example, it has been proposed that the oscillatory behavior of 
the NF-kB transcription factor is key to its capability to integrate cellular signals of 
various origins and to channel them to the appropriate target [SO. 

For this, it is essential to take into account the constraints under which evolu- 
tion has tinkered over the years towards the molecular network studied [7]. Most 
importantly, a biological function not only has to generate an appropriate response 
to a given signal but to do so under predictably or unpredictably varying condi- 
tions, with cross-talk from other functions, while maintaining exquisite sensitivity. 
The combined requirements of robustness and flexibility have undoubtly shaped the 
architecture of cell molecular networks. Accordingly, mathematical modeling has 
been harnessed to gain insight into the dynamical ingredients underlying such prop- 
erties, which make circuits in living cells so different from a random dynamical 
system (see, e.g., ll8TflTT0. 

In this respect, circadian biology is a promising field lfT2T - Tl4"l . Orcadian clocks 
have a well defined function: they keep the time of the day so that the daily changes 
in the environment caused by Earth rotation can be anticipated [fT51[T6ll . These 
biochemical oscillators, made of interacting genes and proteins which interact so 
as to generate oscillations with a period of approximately 24 hours, precisely syn- 
chronize with the day/night cycle so that stable and regular molecular ticks are 
scheduled. Delivering such signals in the fluctuating environment of a cell is a 
formidable challenge in itself, but this task is further complicated by the fact that 
different expression profiles must be generated as the day /night cycle varies across 
the year [T71 - fT9ll , with short days in winter and long days in summer, so as to track 
important moments of the day CUll . Fortunately, while understanding clock ro- 
bustness is a difficult problem, assessing it only requires monitoring the circadian 
oscillator phase. 

To understand the behavior of a functional circadian clock, the autonomous bio- 
chemical oscillator must be considered together with the day /night cycle driving it, 
as it is this very interaction that ensures precise time keeping. The forcing is gen- 
erally parametric: one or several parameters of the internal oscillator are modulated 
by the external cycle, daylight being the principal cue. For example, a clock protein 
may be stabilized by light or degraded faster in the dark. For strong enough cou- 
pling, the phenomenon of entrainment is observed: the period of the oscillator locks 
exactly to that of the external forcing, ensuring that a 24-hour rhythm is generated, 
and a definite phase relationship between the two cycles is maintained H 18112 11 - 1241 . 
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Thus, understanding circadian clock robustness requires identifying the needed 
ingredients both in the autonomous oscillations of the internal oscillator and in its 
response to the external driving cycle. The former question has received much at- 
tention, with special interest in the influence of intrinsic noise (251 - I271 . due to the 
small number of molecules participating in the dynamics, and of temperature vari- 
ations j8l|9l|28l|29l. In contrast, much less effort has been devoted to understand- 
ing how circadian clocks cope with fluctuations in forcing, even though daylight 
strongly fluctuates in natural conditions, throughout the day and from day to day. 
Thus, the same forcing cycle that synchronizes the clock to Earth rotation may 
also reset it erratically. Nevertheless, it has been recognized that daylight fluctu- 
ations affect entrainment accuracy (30l and recently, entrainment to natural L/D 
cycles has been studied more closely (3Tll3~2"l . In particular, it has been hypothe- 
sized that the need to adapt simultaneously to seasonal and weather changes has 
driven clock evolution towards complex architectures featuring several interlocked 
feedback loops Oil . 

Quite unexpectedly, a simple solution to this complex problem was suggested 
by the recent study of the circadian clock of Ostreococcus tauri. This green unicel- 
lular alga displays an extremely simple cellular organization H331I34I and a small and 
compact genome, with low gene redundancy (351 . It shows circadian rhythms in cell 
division [36J, and a genome-wide analysis in L/D cycles revealed rhythmic expres- 
sion for almost all genes, with strong clustering according to biological process (371 . 
Two orthologs of central Arabidopsis clock genes TOC1 and CCA1 were identified 
in Ostreococcus genome by Corellou et al. (381. Overexpression/antisense experi- 
ments supported the hypothesis that TOC1 activates CCA1, which in turn represses 
TOC1, with biochemical evidence of CCA1 binding directly to TOC1 promoter. 

This led Thommen et al. (39l to investigate whether a simple mathematical 
model of a TOC1ICCA1 transcriptional negative feedback loop was consistent with 
microarray data recorded in 12:12 L/D cycles. The agreement between experimen- 
tal data and a simple model based on four differential equations was excellent, pro- 
viding additional support for the TOC1ICCA1 loop hypothesis. Surprisingly, there 
was no signature of a coupling to light in experimental data (391 . as best adjustment 
was obtained with a free-running oscillator model. This finding was confirmed later 
by adjusting simultaneously the microarray data and translational reporter data from 
a different experiment with the same simple mathematical model (401 

Thommen et al. (391 solved this paradox by noting that a precisely timed light 
coupling mechanism can have almost no effect on the oscillator. They exhibited 
several examples of coupling schemes such that an arbitrary parameter modulation 
inside a specific time window induces negligible deviation from the free-running 
profile. Such couplings are invisible when the clock is on time but can reset it 
very efficiently when it drifts, because coupling is then active at a time where the 
oscillator is responsive. 

However, delicate tuning of coupling to light is not only required for shielding 
the clock from weather fluctuations but, perhaps more importantly, also to generate 
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varying clock signals across the year. A natural question is then whether the po- 
tentially conflicting requirements of robustness to fluctuations and of flexibility in 
the temporal profiles generated can be reconciled and satisfied simultaneously. Is 
robustness achieved at the cost of flexibility, or vice versa? For example, simply 
phase- shifting a fixed fluctuation-resistant limit cycle throughout the year would 
allow to track only one moment of the day. 

In this paper, we advance this question by analyzing translational reporter data 
of the Ostreococcus clock genes TOC1 and CCA1 genes recorded for photoperiods 
varying between 2 to 22 hours. Quite remarkably, we observe a complex variation 
of time profiles with photoperiod. However, we also find that all time profiles can be 
accurately adjusted by the simple uncoupled clock model used previously Il39ll40ll . 
provided we assume that the static control parameters of this model can vary with 
photoperiod. Based on our previous work ||39l , this suggests a clock that is robust 
to daylight fluctuations for all photoperiods. 

That such a simple clock, with a clearly identified TOC1-CCA1 one-loop cen- 
tral oscillator, is both robust and flexible is a fundamental and surprising result. It 
suggests that the coupling to light of the TOC1-CCA1 loop is under control of a 
delicately architectured web of yet unidentified additional feedback loops, which 
operate at different scales so as to tune separately the phase dynamics required for 
daily operation, on one hand, and the expression profile modifications required to 
adapt to changing conditions across the year, on the other hand. 

2 Results 

2.1 Construction of the experimental target profiles 

Our analysis is based on time series data of translational reporters of the clock 
genes TOC1 and CCA1, recorded for various light/dark sequences simulating the 
seasonal variations in day and night length (Fig. [[])■ They have been obtained in 
the same conditions as those analyzed by Troein et al. in their modeling study of 
Ostreococcus clock PTTl . We did not use the transcriptional reporter data which are 
also available (Methods). 

Cells were subjected to 24-hour cycles of alternating phases of light (L) and 
darkness (D), termed thereafter day and night. Two experiments were conducted, 
where day lengths for different cell cultures ranged from 2 to 12 hours (Fig.[T]A-F), 
and from 12 to 22 hours (Fig.[T]G-L), respectively. For each cell culture, the lighting 
protocol began with two cycles where photoperiod (day length) was 12 hours (0- 
48 hours), followed by three cycles (termed below photoperiodic cycles) with a 
photoperiod between 2 and 22 hours (48-120 hours), the last cycle being under 
constant illumination (120-144 hours). The two time series with a photoperiod 
of 12 hours (LD 12:12), one in each experiment, have similar profiles although the 
timings differ slightly. Time series corresponding to different cell cultures subjected 
to the same lighting sequence are very consistent, showing the high reproducibility 
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Figure 1: Photoperiodic response of core clock components. Ostreococcus 
TOO: LUC (red) and CCA1:LUC (blue) translational reporter lines were entrained 
under 24-hour day/night cycles with varying day length (A: 12h; B: lOh; C: 8h; D: 
6h; E: 4h; F: 2h; G: 12h; H: 14h; I: 16h; J: 18h; K: 20h; L: 22h). Yellow and black 
bars at the top of each panel indicate periods of light and darkness, respectively. 
Averages of translational reporter lines are drawn as black solid lines. The time 
interval between 96 and 120 hours contains data used for model adjustment. 



of the data, and motivating a quantitative description of experimental data by a 
mathematical model based on biochemical kinetics. 

Our goal is to better understand how Ostreococcus clock entrains to different 
cycles across the year. How do the entrained temporal profiles of the two clock 
genes TOC1 and CCA1, and in particular their peak timings, vary with photope- 
riod? Is the apparent invisibility of coupling previously evidenced in LD 12:12 
experiments [|39ll40l observed for all photoperiods? Therefore, we focus here on 
reproducing data from the third photoperiodic cycle (between 96 and 120 h), as- 
suming that Ostreococcus clock has by then acclimated to the photoperiod change 
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applied at time 48 h, and displays its nominal response to the entrainment cycle. 
This 24-hour interval of time is indicated with a two-head gray arrow in Fig. [TJ 

The time series shown in Fig. \T\ give us the total luminescence emitted by a 
cell population. As such, they may reflect the average single-cell clock dynamics 
as well as variations in the total number of cells or their spatial distribution. For 
example, cell cultures are expected to grow faster in long days since the number 
of cell division events per day was shown to increase with day length [|42l . In 
contrast to this, we expect that the shape of the temporal expression profiles, and in 
particular their peak timings, essentially provides information related to the clock 
dynamics. To avoid capturing amplitude variations possibly biased by population 
dynamics, we normalize the data for the different photoperiods so that they have the 
same maximum expression level. 

In order to keep the mathematical model as simple as possible and avoid over- 
fitting the data, we moreover neglect the fact that luminescence is generated by re- 
porter genes inserted into the genome in addition to the wild ones (possibly inducing 
overexpression effects). We use the translational reporter time series as indicators 
of native protein concentration, an approximation which has been checked in our 
previous works 113911401 and which is acceptable when the photon emission time is 
much larger than the protein degradation time. In fact, a fully detailed model would 
have the simple model described below as limiting case, and would only be needed 
to improve adjustment. Finally, we remove the floor level bias evidenced by Morant 
et al. [|40ll to obtain the target profiles. 

2.2 Orcadian phases of TOC1 and CCA1 display a complex vari- 
ation with photoperiod 

Before carrying out any adjustment, the analysis of the target profiles reveals a com- 
plex orchestration of CCA1 and TOC1 expression by Ostreococcus clock, despite 
its apparent simplicity. This is illustrated by Fig. [2l which displays the time posi- 
tions of TOC1 and CCA1 expression peaks as a function of Zeitgeber Time (ZT), 
measured from last dawn. It appears clearly that the TOC1-CCA1 oscillator does 
not respond to varying photoperiod by simply globally shifting its phase but relies 
on a mechanism that controls differentially the peak timings of CCA1 and TOC1. 
TOC1 peak tends to track dusk except for very short days where its timing becomes 
roughly fixed relative to dawn. For long days, CCA1 peak tracks dawn, occurring 
around ZT21, while for short days it combines information about both dawn and 
dusk, occurring approximately one hour after the middle of the night. 

The time delay between TOC1 and CCA1 expression peaks, as well as their 
times of occurrence in the lighting cycle, are therefore critical clock features. One 
objective of this work is to check whether the variations of TOC1 and CCA1 tem- 
poral profiles with photoperiod, including interpeak delay and peak widths, can be 
reproduced by a mathematical model. Interestingly, it is for the pivotal case of LD 
12:12 that the TOC1-CCA1 delay is larger and the expression peaks narrower. As 
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the interpeak delay decreases, above or below a photoperiod of 12 hours, expression 
peaks tend to broaden. 

For each photoperiod, the relative position of the two expression peaks remains 
compatible with the hypothesis of a two-gene loop where TOC1 activates CCA1 
and CCA1 represses TOC1, except for a photoperiod of 22 hours, where the two 
peaks coincide. Our assumption that luminescence signals reflect true protein con- 
centrations reaches then its limits. Recall that luminescence data of CCA 1 -LUC and 
TOC1-LUC come from different cell lines. Insertion of CCA1-LUC (resp. TOC1- 
LUC) induces an overexpression of CCA1 (resp. TOC1), which can be shown to 
shorten (resp. lengthen) the free-running period (FRP). These opposite FRP varia- 
tions result in slight antagonist phase shifts of expression profiles, leading to over- 
lap of the expression peaks. This problem could be worked around simply by using 
more detailed models of the two transgenic clocks (including the CCA1-LUC and 
TOC1-LUC genes), however we will see later that this would only be needed for 
extreme photoperiods, encouraging us to preserve the simplicity of our model. 



CM 




4 8 12 16 20 24 

Zeitgeber time (h) 

Figure 2: Complex response of Ostreococcus clock to photoperiod changes. Tim- 
ing (ZT) of the concentration peaks of TOCl:LUC (red) and CCA1:LUC (blue) are 
shown as a function of day length. The yellow and gray backgrounds indicate inter- 
vals of light and darkness, respectively. Dashed lines correspond to middle of the 
day and middle of the night. 
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2.3 Mathematical model 



The minimal model of the TOC1-CCA1 transcriptional feedback loop consists of 
the following four ordinary differential equations: 

»> A,j K M jM T 

M T = u T -i ; Omt 1 ( ia ) 

l + (P c /Pco) nc t K Mt +M t 

Kp Pt 

P t = p T M T -8p T tr ^ l D (lb) 
&p T + Pt 

,> , MPt/Pto) ,1t ~ K Mc M c 

Mr = A/rH ~, — —, s OMr (lc) 

PC l + (P T /P T o) nT Mc K Mc +M c 

P C = PcM c -5p c -^^- (Id) 
K P C + "c 

Eqs. (OQ) describe the time evolution of mRNA concentrations Mq and Mj and pro- 
tein concentrations Pq and Pj for the CCA1 and TOCi genes, respectively, as they 
result from mRNA synthesis regulated by the other protein, translation and enzy- 
matic degradation. TOC1 transcription rate varies between /jj at infinite CCA1 
concentration and /uj + Xj at zero CCA1 concentration according to the usual gene 
regulation function with threshold Pqq and cooperativity hq- Similarly, CCA1 tran- 
scription rate is /uc (resp., /uc + ^c) a t zero (resp., infinite) TOC1 concentration, 
with threshold Pjo and cooperativity nj. TOC1 and CCA1 translation rates are 

and (3c, respectively. For each species X, the Michaelis-Menten degradation 
term is written so that dx is the low-concentration degradation rate and Kx is the 
saturation threshold. Their expression is mathematically equivalent to the usual for- 
mulation using the Michaelis-Menten constant Kx and the maximum degradation 
speed V max = d x K x . 

We had found previously that Eqs. ® with no parametric modulation, corre- 
sponding to a free-running oscillator, could reproduce simultaneously microarray 
and translational reporter data recorded under two different LD 12:12 experiments, 
with excellent agreement Il40l . To evaluate adjustment reproducibility, we checked 
whether Eqs. (OQ) using the best-fitting parameter set obtained in HOl could repro- 
duce the two LD 12:12 datasets used here. As Fig. [3] shows, there is an excellent 
agreement with the short-day experiment recording, which is the third experiment 
adjusted by this model with this parameter set! For the long day experiment, there 
is a noticeable shift in the CCA1 peak, which is not understood. As we see below, 
however, this will not affect our main conclusions. 



2.4 Photoperiod-dependent free-running oscillator models ad- 
just experimental data precisely 

A simple way to obtain a robust coupling for all photoperiods would be to globally 
phase shift the robust oscillator evidenced in previous work 1139114011 depending on 
photoperiod, keeping the two profiles and their separation unchanged. However, the 
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very fact that T0C1-CCA1 delay and peak widths vary significantly with photope- 
riod totally excludes such possibility. It is then quite mysterious how robustness can 
be achieved simultaneously for each of the different profiles observed. 

In order to gain insight into how robustness combines with the flexibility ob- 
served, we tested adjustment of the expression profiles profiles by free-running 
TOC1-CCA1 oscillator models whose parameters are allowed to vary with pho- 
toperiod. This amounts to carrying out for each photoperiod the same analysis that 
we applied to a LD 12:12 dataset in our previous works [|39ll40ll . Parameter values 
of modelQ]are optimized to adjust this model to the experimental profiles under the 
constraint that the FRP is 24 hours (as discussed in [|39l . this technical simplification 
does not imply that the actual FRP is exactly 24 hours). 

Fig|4]-A shows clearly that an impressive agreement between experimental and 
numerical profiles can be obtained for all photoperiods. In particular, the timings 
of the TOC1 and CCA1 expression peaks are very well reproduced (Fig.|4]-B). By 
applying the same procedure to random target profiles, we found that the probability 
of obtaining such a good agreement for all photoperiods by chance is exceedingly 
low, well below 10 -4 (Methods). Thus, this finding strongly suggests that a strong 
evolutionary constraint has shaped the expression profiles of Ostreococcus clock. If 
our interpretation of the invisible coupling behavior is correct, this constraint is the 
necessity of maintaining robustness to daylight fluctuations all across the year. 

That Ostreococcus clock generates time profiles which remain so close to that of 
a free-running oscillator while strongly varying with photoperiod is a very impor- 
tant result. It shows that robustness and flexibility can be simultaneously achieved 
in a simple clock built around a simple one-loop, two-gene, core oscillator. This 
flexibility appears clearly in Fig.|H-C, which shows the projections into the TOC1- 
CCA1 plane of the different limit cycles observed for all photoperiods. 

For each photoperiod, excellent adjustment was obtained in a wide region of pa- 




Zeitgeber time (h) 

Figure 3: Comparison of data recorded in LD 12:12 with numerical solutions of 
model ([T]) using the parameter set obtained in Ref. Il40ll . TOC1 (resp., CCA1) pro- 
files are drawn in red (resp., blue). (A) Short-day experiment, corresponding to 
signal in Fig.l-A; (B) Long-day experiment, corresponding to signal in Fig.l-G. 
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rameter space (see [|39ll for more details), a phenomenon which is often observed in 
systems biology models Il43l . Independent adjustments for each photoperiod there- 
fore yield wildly varying parameter values as a result. We therefore used a modified 
goodness of fit estimator penalizing parameter value dispersion across photoperi- 
ods, to ensure that parameter value variations across the photoperiod range were 
only as large as needed to reproduce the different time profiles. Fig.[5Jshows the ob- 
tained parameter values vary as a function of photoperiod. For most parameters, a 
limited and smooth variation with photoperiod is observed, consistent with the idea 
of a tuning of the TOC1-CCA1 loop depending on photoperiod. The parameter val- 
ues for the two LD 12:12 experiments are generally relatively consistent between 
themselves as well as with the parameter values obtained in ll39l . 

It is difficult to assess the biological significance of these curves at this stage 
because time series are purposely adjusted up to a scaling both in concentrations 
and in time (to keep the FRP at 24 hours). While this guarantees that only the clock 
dynamics is probed, possible systematic variations in expression levels with pho- 
toperiod are not captured and should manifest themselves as systematic variations 
in the best fitting parameter values (Methods). It is thus plausible that less parame- 
ters do vary with photoperiod than appears in Fig. [51 However, we concentrate here 
on our strategy to evidence the free-running behavior of the TOC1-CCA1 oscillator 
at all photoperiods, leaving precise parameter estimation for a future work. 

2.5 Robustness to daylight fluctuations 

Here we illustrate how robustness to daylight fluctuations can be achieved with a 
light coupling mechanism that does not visibly affect the dynamics of the core os- 
cillator in entrainment conditions. As an example, we will assume that the core 
TOC1-CCA1 oscillator is driven by the day/night cycle via a transient modulation 
of parameter dp T (TOC1 degradation rate). For the sake of simplicity, this modu- 

Table 1 : Parameter values obtained by adjusting experimental data with a free- 



running model with 16 photoperiod-dependent parameters 





Ref 40 


Short day experiments 




Long day experiments 


Day length 


12 


2 


4 


6 


8 


10 


12 




12 


14 


16 


18 


20 


22 


RMS Error 


7.85 


7.3 


4.5 


3.7 


3.5 


3.8 


4.5 




5.3 


4.6 


5.9 


7.4 


9.0 


9.7 


iif 


2 


2 


2 


2 


2 


2 


2 




2 


2 


2 


2 


2 


2 


fi c (nM.rr 1 ) 


0.153 


0.262 


0.174 


0.216 


0.163 


0.227 


0.197 




0.626 


0.379 


0.787 


0.871 


0.622 


0.235 


Xc (nMrT 1 ) 


3.10 


7.43 


3.56 


3.49 


2.54 


2.98 


2.45 




3.04 


3.34 


3.99 


3.56 


2.72 


2.46 


Pro (nM) 


18.7 


48.1 


38.8 


38.5 


35.4 


36.8 


26.8 




37.0 


42.6 


26.3 


26.6 


21.8 


16.2 


Pc(h _1 ) 


2.83 


4.90 


4.52 


4.4X 


5.71 


5.34 


4.77 




5.72 


10.3 


7.88 


9.29 


13.1 


15.5 


"c 


2 


2 


2 


2 


2 


2 


2 




2 


2 


2 


2 


2 


2 


fij (nM.h~ ! ) 


0.467 


2.54 


2.82 


2.50 


2.58 


2.59 


2.96 




2.54 


2.54 


3.20 


3.20 


3.22 


2.92 


X T (nM.ri" 1 ) 


487 


108 


86.0 


84.5 


108 


84.9 


86.0 




101 


79.8 


93.9 


93.2 


90.2 


94.1 


Pea (nM) 


4.51 


14.9 


8.06 


7.90 


6.52 


5.57 


7.24 




4.89 


6.78 


5.11 


4.87 


4.91) 


6.00 


p> (h- 1 ) 
1/8m c (h) 


0.812 


0.072 


0.105 


0.135 


0.185 


0.457 


0.801 




0.157 


0.226 


0.086 


0.063 


0.037 


0.024 


0.195 


0.137 


0.193 


0.205 


0.232 


0.183 


0.209 




0.245 


0.274 


0.245 


0.275 


0.423 


0.376 


l/5 Pc (h) 


2.36 


2.91 


1.86 


1.68 


1.68 


1.83 


2.08 




1.16 


0.654 


0.444 


0.370 


0.286 


1.79 


1/&M T (h) 


0.129 


0.074 


0.097 


0.098 


0.103 


0.114 


0.119 




0.106 


0.107 


0.115 


0.120 


0.125 


0.122 


1/Spj. (h) 


0.199 


0.538 


0.644 


0.630 


0.400 


0.318 


0.273 




0.394 


0.339 


0.533 


0.563 


0.479 


0.105 


K Mc (nM) 


0.407 


0.303 


0.250 


0.250 


0.268 


0.274 


0.273 




0.345 


0.303 


0.407 


0.463 


0.559 


3.364 


Kp c (nM) 


75.9 


2237 


1386 


1258 


1168 


1105 


843 




1579 


2239 


12054 


13648 


14678 


10575 


*M T (nM) 


28.3 


8.25 


8.17 


8.09 


9.21 


6.70 


7.72 




6.68 


6.65 


7.22 


6.93 


7.51 


5.70 


Kp T (nM) 


2.76 


2.50 


2.98 


3.04 


2.99 


3.71 


3.31 




2.91 


2.79 


3.24 


3.12 


1.70 


0.273 
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Figure 4: Adjustment of TOC1 and CCA1 profiles by a free-running oscillator 
model. (A) Measured data shown in Figure 1 are compared with numerical so- 
lutions of a free-running oscillator model (Eqs. (OQ)) where all parameters have been 
adjusted for each photoperiod. (B) Comparison of the ZT peak timings of TOC1 
(red circles) and CCA1 (blue squares) expression, for experimental data (full sym- 
bols) and model (empty symbols). An excellent agreement is observed both for 
profile shapes and timings. (C) Projection of the limit cycle into the TOC1-CCA1 
plane for the different photoperiods. Line color indicates photoperiod, and ranges 
from dark blue (2-hour day length) to dark red (22-hour day length), light blue and 
light red corresponding to short-day and long-day 12:12 protocols, respectively. 



lation occurs in a temporal window which is fixed relative to the day /night cycle. 
There is no loss of generality in doing so if we only intend to understand the behav- 
ior in entrainment, not the dynamics of resetting after a large phase excursion. How 
such a parametric forcing ensures a stable phase relation between the clock and the 
day/night cycle can be understood as follows. 

The action of a parametric modulation on an oscillator is described by the im- 
pulse phase response curve (iPRC) [8l[32]|44]|, which gives the phase shift induced 
by a short parameter perturbation with unit time integral as a function of the position 
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Figure 5: Values of the best-fitting parameters used in Fig. |4]as a function of pho- 
toperiod. The red circle indicates the parameter value adjusted in Ref. [40J. 



along the limit cycle. Fig. [6]- A shows the iPRCs of TOC1 degradatation modulation 
computed for the entrained limit cycles observed at varying photoperiods. Although 
the shape of these iPRC slightly varies for these distinct limit cycles, it is remark- 
able that the fiat part of the iPRC persists. This conserved feature of the iPRC, also 
known as a dead zone, reflects an insensitivity of the limit cycle to TOC1 degrada- 
tion modulation at the begining of the day. This implies that the coupling profile 
can be designed so as to match these intervals of insentivity, endowing the clock 
with robust entrainment by fluctuating forcing [44J. 

It is to note that, in the limit of weak coupling, the total phase shift 8(|) induced 
by parameter perturbation profile 8p(t) is given by 8(|) = f t off Z(t)8p(t)dt where Z(t) 
is the iPRC and t on and f ff are the start and end times of the coupling window. For 
a FRP of approximately 24 hours, stationary operation requires 8(|) ~ so that the 
coupling window [? n,^off] must be located around a zero of the iPRC Z{t) (a stable 
operation also requires a negative derivate of Z[t) at this zero). 

For each photoperiod, we searched for coupling windows inside which TOC1 
degradation can be increased without degrading model adjusment compared to the 
uncoupled model. This is in principle a stringer requirement than the zero phase 
shift condition since this also implies that the deviation from the limit cycle is min- 
imal. Fig.[6]-B shows the locations of such windows where coupling activation was 
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Figure 6: Robust coupling to light at all photoperiods. (A) Phase Response Curves 
(iPRCs) of the different photoperiod-dependent free-running oscillators in response 
to an impulse perturbation of TOC1 degradation rate. (B) Coupling windows in 
which TOC1 degradation rate is multiplied by two and beginning at ZT0 such that 
adjustment of model to experimental data is preserved (no deformation of the limit 
cycle). (C) Entrained TOC1/CCA1 oscillatory profiles are insensitive to fluctuations 
of TOC1 degradation rates (modulation factor uniformly distributed between 1 and 
3) inside the coupling windows in Panel B. (D) Limit cycles corresponding to the 
profiles shown in Panel C, showing the small residual fluctuations. 



fixed at ZT0 and increase of TOC1 degradation rate during the coupling window at 
100 %. In fact, the timings of these windows almost did not depend on the mod- 
ulation factor chosen, indicating that the same limit cycle, close to the uncoupled 
one, is obtained at different coupling strengths (hence light levels). Interestingly, 
window timings are relatively symmetrical around the LD 12:12 protocol. 

As Fig. [6J-C shows, this insensitivity of the adjustment with respect to mod- 
ulation strength ensures that the circadian oscillator is robust to daylight fiuctua- 
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tions. We subjected the limit cycles obtained for three different photoperiods (LD 
12: 12, LD 6: 18 and LD 18:6) to a random sequence of modulation depths. On each 
day, TOC1 degradation rate was multiplied inside the coupling window shown in 
Fig.[6]-B by a random factor uniformly distributed between 1 and 3. This describes 
adequately the effect of daylight variations from one day to the next, which are the 
most disruptive because they are resonant with the oscillator 11391 . It can be seen that 
the circadian clock ticks very robustly in the three cases, delivering signals almost 
identical to those of a free-running oscillator even though there is a strong coupling 
operating for several hours (Fig.[6]-B). Obviously, this only holds when the oscillator 
has the phase observed experimentally. Otherwise, the zero phase shift condition is 
not satisfied and the coupling rapidly resets the clock to the correct time. 

Fig. |6]-D displays the three limit cycles under random forcing, showing the small 
residual fluctuations. Interestingly, these are stronger for the two extreme photoperi- 
ods than for the LD 12:12 protocol which is impressively insensitive to fluctuations. 
This figure, with its clearly distinct limit cycles, illustrates nicely how robustness 
to daylight fluctuations and flexibility can be combined in a simple clock provided 
coupling to light is suitably designed. 

To summarize this section, adjustment of circadian signals by a free-running 
oscillator model does not imply that there is no coupling, but rather that coupling is 
scheduled so that it does not affect the oscillator when the clock is entrained, thereby 
shielding it from daylight fluctuations. Importantly, it was shown in 11391 that even 
when the free-running period is different from 24 hours, a similar mechanism is 
possible. Coupling then acts to rescale the oscillation period to 24 hours, but leaves 
no other print besides a small amplitude variation, as the waveforms can be adjusted 
by a free-running oscillator model. 

3 Discussion and Conclusion 

In this work, we have studied the response of Ostreococcus circadian clock to 
day/night cycles of various day lengths, such as those induced by seasonal changes. 
The analysis of the photoperiod-dependent expression profiles of the two central 
clock proteins TOC1 and CCA1 allowed us to unveil two remarkable features of 
Ostreococcus clock photoperiodic behavior. 

First, a flexible and complex variation of expression profiles with photoperiod is 
observed, reminiscent of that observed in the higher plant Arabidopsis 11201 . There is 
not only a complex response of the oscillator phase to varying photoperiod [1411 but 
also the time interval between the two expression peaks and the waveform of the 
profiles change significantly with photoperiod. More precisely, TOC1 expression 
tends to follow the day-night transition, except in shorter days where it remains 
fixed relatively to dawn, whereas CCA1 expression is scheduled in the middle of the 
night, except for longer days where it tracks dawn (Fig. |2]). Therefore, at least two 
phases rather than a single one are needed to describe entrainment of Ostreococcus 
clock to different photoperiods. More generally, it is consistent with the view that a 
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circadian clock must deliver different signals in different seasonal contexts because 
different physiological processes must be controlled at different times of the day, 
which is indeed the case in Ostreococcus 113711 . 

Second, we found that protein expression data for all photoperiods can be ad- 
justed with surprising accuracy by a minimal model of the TOC1-CCA1 transcrip- 
tional loop with no driving by light, provided that its kinetic parameters are allowed 
to vary with photoperiod to account for the flexible response observed. This re- 
markable finding not only confirms the central role of the TOC1-CCA1 loop in 
Ostreococcus circadian clock ll39l - l4Tll45l but also indicates that the light coupling 
mechanisms which synchronize the clock become invisible when the clock is en- 
trained. As we showed, a simple explanation for this phenomenon is that activation 
of the light input pathway coincides with a "dead zone" of the phase response curve 
of the circadian oscillator Il44l . Such a design ensures that the same profiles are 
generated at different daylight intensities, in which case random fluctuations have 
no effet. This finding suggests that molecular signals from this clock have been 
strongly shaped by the requirement of being robust to daylight fluctuations all along 
the year. This constraint is all the stronger in Ostreococcus as the light perceived 
by this marine organism varies not only due to sky cover but also depending on 
distance to surface and water turbidity. 

Flexibility and robust synchronization to the day/night cycle both rely on the 
architecture of the light input pathways and feedback loops interacting with the 
TOC1-CCA1 oscillator, about which little is currently known despite recent progress Il46l 
1471 . Together, our results suggest that there are there is a clear distinction of light 
coupling processes according to dynamical role and time scale. 

On the one hand, gated coupling mechanisms, where a kinetic parameter is mod- 
ulated for a few hours, are required to ensure entrainment and keep the oscillation 
phase stable from day do day. Gated light inputs is a feature of many circadian 
clocks and it has proposed that they are critical for generating appropriate timings 
under different photoperiods [18]. Closely related phenomena are light adaptation 
or response saturation H32[|48[|49l . We have proposed that in a robust clock, the fast 
synchronizing inputs must be tuned so that they do not deform the entrained limit 
cycle, which would otherwise be exposed to daylight fluctuations [|39l . Therefore, 
the flexible variation of expression profiles with photoperiod needed to adapt clock 
signals to different seasons must entirely rely on slow inputs modifying kinetic pa- 
rameters of the core oscillator on a time scale of several days. Indeed the existence 
of such slow kinetic changes naturally explains why the variation of expression pro- 
files with photoperiod is perfectly matched by a free-running oscillator model with 
photoperiod-dependent parameters. Such slow changes may be related to light accu- 
mulation processes, coincidence mechanisms or metabolism, and more generally to 
the presence of molecular actors constitutively expressed at photoperiod-dependent 
levels. In fact, this separation in fast mechanisms maintaining phase and slow mech- 
anisms controlling expression profiles is only a necessary condition for robustness, 
as the synchronizing couplings and their phase response curves must also satisfy 
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specific contraints to ensure a reproducible phase dynamics over a wide range of 
operating conditions Il44l . 

Admittedly, implementing the design principles we have unveiled in a detailed 
mechanistic model of Ostreococcus clock is not an easy task, if only because of the 
number of slow inputs needed to reproduce the variation of expression profiles with 
photoperiod. In this respect, it would be interesting to analyze the Ostreococcus 
clock model proposed by Troein et al. IflTl in the light of our results. This model, 
which takes into account the detailed dynamics of the luminescent reporters, was 
constructed from an extensive dataset obtained in essentially the same experimen- 
tal conditions as our time series. Both limit cycle profiles and transient behavior 
(change in photoperiod, transition into constant light) were taken into account to 
adjust the mathematical model. Troein et al. concluded that the complex behav- 
ior of the one-loop Ostreococcus clock could be explained by the presence of five 
different light inputs into CCA1-TOC1 loop, seemingly contradicting our results. 
It may be that these light inputs are gated by the profiles of the actors they affect. 
Interestingly, the profiles generated by this model are relatively smooth, except at 
dawn and dusk, which could indicate a weak effective coupling. Otherwise, the 
light accumulator present in this model is similar to a slow input, although it seems 
to be operating on a time scale too short to ensure robustness, and probably cannot 
reproduce alone the flexibility observed in expression profiles. 

It is also interesting to compare our results to Ref. EOl where the flexibility of 
increasingly complex model of Arabidopsis clock was studied and compared to ex- 
perimental results. It was shown that mechanistic models of this circadian require at 
least three feedback loops to generate timing patterns that are not limited to tracking 
dawn or dusk and which reproduce those observed experimentally. This finding is 
consistent with the theoretical prediction that flexibility requires several interlocked 
feedback loops j8l|9l[5Ql|. Yet, it seems that a one-loop model with photoperiod- 
dependent parameters is able to reproduce very complex patterns where, for exam- 
ple, CCA1 can track dawn or the middle of the night depending on photoperiod. 
This suggests that more attention should be given to the role of slow light inputs, 
serving as photoperiod sensors. Such slow effects of light have been recently evi- 
denced in mice 11511 . It is probably the case that in Ostreococcus, these slow inputs 
are controlled by additional feedback loops that remain to be discovered. What 
complicates their identification is that they seem to be invisible in entrainment con- 
ditions by design. 

Clearly, Ostreococcus clock has not yet revealed all its secrets. Our partial re- 
sults suggest that unraveling them may unveil important design principles in circa- 
dian biology, thanks to the surprising agreemens between mathematical models and 
experimental data that can be obtained in this model organism. Whatever molecular 
model of Ostreococcus clock eventually emerges, an important lesson from Ostreo- 
coccus clock response to photoperiod changes will undoubtly be that even a simple 
circadian clock can combine robust and flexible mechanisms to adapt to changing 
weather and seasons. 
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4 Methods 



4.1 Experimental data 

The experimental data were obtained as described by Corellou et al. Il38ll . in essen- 
tially the same experimental conditions as used by Troein et al. [HTTl . We did not use 
the transcriptional reporter data available. Indeed, Djouani-Tahri et al. ll46Tl showed 
directly, using inhibitors of translation, that the time between luciferase synthesis 
and photon emission is extremely long, of the order of 8-10 hours. This is the 
time required to observe a significant decrease in luminescence when translation is 
blocked. A conclusion of our previous works H39U40H is that transcriptional activities 
of both TOC1 and CCA1 are confined to much shorter time interval. As a result, the 
signal of interest, namely transcriptional activity, is probably averaged out in the 
transcriptional reporter luminescence data, which may destabilize the adjustment 
procedure 

4.2 Target profiles 

The target profiles of protein concentration are obtained from the luminescence 
records of the translational reporter lines during the third cycle of photoperiod (from 
96 h to 120 h in FigQ])- Luminescence time series display variations in amplitude 
from day to day due to fluctuations in the number of cells contributing to light 
emission and other unknown factors. To correct this effect approximately and obtain 
periodic target profiles, the time series were divided by a first order polynomial 
function of time interpolating between the luminescence intensities at 96 h and 120 
h. The floor level was then removed from the time series to correct for the bias 
evidenced in [|40] and periodic profiles for all photoperiods were then rescaled to 
have the same maximum value. 

4.3 Adjustment and goodness of fit 

Model (QQ) has 16 free continuously varying parameters besides the cooperativities 
nc and nj which are set to the integer values 2. Two solutions of Eqs. (OQ) that 
have the same waveforms up to a scale factor are considered equivalent. Although 
this allows us to factor out population effects, this makes it more complicated to 
compare parameter sets resulting from optimization, because each one is actually 
the representative of a four-dimensional manifold of equivalent points in parameter 
space. This is adequate for our purposes since our goal here is to establish the 
relevance of the free-running TOC1-CCA1 oscillator model for all photoperiods, 
not to estimate accurate parameter values. To measure the goodness of fit for a 
given parameter set, the two numerical protein profiles are rescaled to have the 
same maximum value as the experimental profiles and compared to the latter by 
computing the root mean square (RMS) error. 
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4.4 Optimization 

Adjustment was carried out by using a population-based metaheuristic method (har- 
mony search [52]) for the initial large-scale search, followed by a nonlinear opti- 
mization procedure based on a Modified Levenberg-Marquardt algorithm (routine 
LMDIF of the MINPACK software suite [53]) to refine the optimal parameter val- 
ues. The procedure constrains the FRP at a value of 24 hours and the goodness of 
fit includes both RMS error and a penalty proportionnal to the Eulerian distance 
between current parameter set and the reference one obtained in [40J. This penalty 
maximizes correlation between the different best-fitting parameter sets and makes 
their comparison easier. Numerical integration of ordinary differential equations 
was performed with the SEULEX algorithm [|54l . The exhaustivity of the harmony 
search stage and the convergence of the adjustment were monitored by checking 
that the optimum was reached repeatedly. 

4.5 Probability of adjustment by free-running oscillators 

In order to show that the simultaneous adjustment of photoperiodic data by free- 
running oscillators is biologically significant, it is important to show that such a 
numerical result cannot be obtained by chance. We therefore generated a large 
number of random target profiles which were then fed to the optimization procedure. 
The random protein profiles (one for TOC1 and one for CCA1) featured a single 
peak per period, obtained by interpolating three control points with a cubic spline. 
A TOC1-CCA1 delay was then chosen at random between 2 and 12 hours. This 
yields smoothly varying profiles similar to the experimental ones and to those which 
can be generated by a free-running 4-ODE model. 

We found that for this set of random profiles, the probability of obtaining a 
RMS error as good as obtained in Fig. [4] is always lower than 0.4. The probability 
of obtaining 11 such adjustments is thus bounded above by 4 10~ 5 . 
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